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Abstract 

The paper collects the various pieces of information concern- 
ing the relative size of m u , and m s . A coherent picture re- 
sults, which constrains the mass ratios to a rather narrow range: 
m u /m d = 0.553 ± 0.043, m s /m d = 18.9 ± 0.8. 
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1. I wish to show that recent results of chiral perturbation theory al- 
low a rather accurate determination of the relative size of m u , rrid and m s . 
The paper amounts to an update of earlier work JTJ-0 based on the same 
method. Sum rules and numerical simulations of QCD on a lattice represent 
alternative approaches with a broader scope - they permit a determination 
not only of the ratios m u : rrid : m a , but also of the individual quark masses, 
including the heavy ones. The sum rule results for the ratios are subject to 
comparatively large errors |7|, §. Concerning the lattice technique, consid- 
erable progress has been made ||, Ho]. It is difficult, however, to properly 



account for the vacuum fluctuations generated by quarks with small masses. 
Further progress with light dynamical fermions is required before the num- 
bers obtained for m u /md or m s /md can be taken at face value. 

2. The quark masses depend on the renormalization scheme. Chiral per- 
turbation theory treats the mass term of the light quarks, m u uu + rriddd + 



m s ss, as a perturbation [JTT], |12| . It exploits the fact that, for mass inde- 
pendent renormalization schemes, the operators uu, dd and ss transform as 
members of the representation (3, 3*) + (3*, 3). Since all other operators with 
this transformation property are of higher dimension, the normalization con- 
ventions for m u , irid and m s then differ only by a flavour- independent factor. 
The factor, in particular, also depends on the renormalization scale, but in 
the ratios m u /m d and m s /m d , it drops out - these represent convention- 
independent pure numbers. 

3. The leading order mass formulae for the Goldstone bosons follow from 
the relation of Gell-Mann, Oakes and Renner. Disregarding the electromag- 
netic interaction, they read M%+ = {m u + md)B, M\+ = {m u + m s )B, 
Mft = (rrid + m s )B, where the constant of proportionality is determined by 
the quark condensate: B = |(0|mu|0)|/-F^. Solving for the quark masses and 
forming ratios, this constant drops out, so that m u /rrid and m s /md may be 
expressed in terms of ratios of meson masses. Current algebra also shows 
that the mass difference between the 7r + and the 7r° is almost exclusively due 
to the electromagnetic interaction - the contribution generated by rrid 7^ m u 
is of order (rrid — ^u) 2 and therefore tiny. Moreover, the Dashen theorem 
states that, in the chiral limit, the electromagnetic contributions to M|- + 
and to are the same, while the self energies of K° and tt° vanish. Us- 
ing these relations to correct for the electromagnetic self energies, the above 
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lowest order mass formulae yield ffl] 



m u _ M| H -M| ■ 2\ri,-M^_ () .. _ 



ms M|o + M| + - M* + 



20.1 



4. These relations and the analysis given below are based on the hypo- 
thesis that the quark condensate is the leading order parameter of the sponta- 
neously broken symmetry. This hypothesis is questioned in refs. [13|, where 



a more general scenario is described, referred to as generalized chiral per- 
turbation theory: Stern and co-workers investigate the possibility that the 
correction of 0(m 2 ) in the expansion M 2 = {m u +m l i)B J r O{m 2 ) is comparable 
with or even larger than the term that originates in the quark condensate. 
Indeed, the available evidence does not exclude this possibility, but a beau- 



tiful experimental proposal has been made ||14|| : 7r + 7r~ atoms decay into a 
pair of neutral pions, through the strong transition 7r + 7r~ — > 7r°7r°. Because 
the momentum transfer nearly vanishes, the decay rate is determined by the 
combination a$ — 0,2 of S-wave irn scattering lengths. Since chiral symmetry 
implies that Goldstone bosons of zero energy do not interact, ao, 0,2 vanish in 
the limit m K ,m,i->0. The transition amplitude, therefore, directly measures 
the symmetry breaking generated by m u , m^. Standard chiral perturbation 



theory yields very sharp predictions for ao, a>2 @, while the generalized sce- 
nario does not |nj. A measurement of the lifetime of a 7r + 7r~ atom would 
thus allow us to decide whether or not the quark condensate represents the 
leading order parameter. 

5. The contributions of first non-leading order were worked out in ref. 
[0. As is turns out, the correction in the mass ratio (M^o — M^+)/(M|- — 
Ml) is the same as the one in M^/M^: 

m k tti ~\~ ttl s , 2 



Ml m u + m d 



{l + A M + 0(m 2 )} , (2) 



M 



K" ~ M K+ _ m d -m u _LrVrr. 2 U 

Ml-Ml -^^T {1 + AM + ° (m)} ' 
The quantity A M accounts for the breaking of SU(3) and is related to the 
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effective coupling constants L5 and L§: 

8(Mj- - 



Mi 



(2L S - L 5 ) + xlogs . (3) 



The term xl°gs stands for the logarithms characteristic of chiral perturbation 
theory (for an explicit expression, see ||12|| ). The above relations imply that 
the double ratio 

_ 2 m 2 -m 2 . s 

is given by a ratio of meson masses, up to corrections of second order, 

Ml Ml -Ml 
Q ~ Ml Ml -Ml + {1 + ° {mn ■ (5) 



The result may be visualized by plotting m s /md versus m u /md ]T7|]. The 
constraint then takes the form of an ellipse, 

m,,\ 2 1 f m H \ 2 , . 

m d J Q 2 \m d J 

with Q as major semi-axis, the minor one being equal to 1 (for simplicity, I 
have discarded the term m 2 /m 2 , which is numerically very small). 

6. The meson masses occurring in the double ratio (|5p refer to pure QCD. 
Correcting for the electromagnetic self energies with the Dashen theorem, the 
quantity Q becomes 

2 (Mj-o + M 2 K+ - Ml + + Ml )(M 2 K0 + M 2 K+ - M 2 n+ - M^ 2 ) 

4Ml Q (M 2 K0 -M 2 K+ + Ml + -Ml ) ■ {l) 

Numerically, this yields Qe> = 24.2. For this value of the semi-axis, the 
ellipse passes through the point specified by Weinberg's mass ratios, which 
correspond to A M = 0, Q = Qd- 

The Dashen theorem is subject to corrections from higher order terms in 
the chiral expansion, which are analysed in several recent papers. Donoghue, 



Holstein and Wyler |18| estimate the contributions arising from vector meson 
exchange and conclude that these give rise to large corrections, increasing the 
value (M K +—M K o) e m ^ = 1.3 MeV predicted by Dashen to 2.3 MeV. According 



to Baur and Urech |B|, however, the model used is in conflict with chiral 
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symmetry: although the perturbations due to vector meson exchange are 
enhanced by a relatively small energy denominator, chiral symmetry prevents 
them from being large. In view of this, it is puzzling that Bijnens [EUJ, who 



evaluates the self energies within the model of Bardeen et al., finds an even 
larger effect, (M K + — M K o) ejn . ~ 2.6 MeV. The implications of the above 
estimates for the value of Q are illustrated on the r.h.s. of fig. 1. 

Recently, the electromagnetic self energies have been analysed within lat- 
tice QCD (T|. The result of this calculation, (M K + - M K o) e , m , = 1.9 MeV, 



indicates that the corrections to the Dashen theorem are indeed substantial, 



although not quite as large as found in refs. |L8], p0| . The uncertainties of the 
lattice result are of the same type as those occurring in direct determinations 
of the quark masses with this method. The mass difference between K + and 
K°, however, is predominantly due to m^>m n , not to the e.m. interaction. 
An error of 20% in the self energy affects the value of Q by only about 3%. 
The terms neglected when evaluating Q 2 with the meson masses are of order 
(M| - M2) 2 /M 4 , where M is the mass scale relevant for the exchange of 
scalar or pseudoscalar states, M ~M ao ~M V ' j4|. The corresponding error 
in the result for Q is also of the order of 3% - the uncertainties in the value 
Q = 22.8 that follows from the lattice result are significantly smaller than 
those obtained for the quark masses with the same method. 

7. The isospin- violating decay i] — > 3tt allows one to measure the semi-axis 
in an entirely independent manner f21 |. The transition amplitude is much less 



sensitive to the uncertainties associated with the electromagnetic interaction 
than the K°-K + mass difference: the e.m. contribution is suppressed by chi- 



ral symmetry and is negligibly small [22]. The decay 77 — > 3tt thus represents 



a sensitive probe of the symmetry breaking generated by rrid — m u . It is con- 
venient to write the decay rate in the form r^_^. 7r + 7r - 7r o = r (Qd/Q) 4 , where 
Qd is specified in eq. (|7|). As shown in ref. |21|, chiral perturbation theory 
to one loop yields a parameter-free prediction for the constant T . Updating 
the value of F n , the numerical result reads T = 168 ± 50 eV. Although the 
calculation includes all corrections of first non-leading order, the error bar is 
large. The problem originates in the final state interaction, which strongly 
amplifies the transition probability in part of the Dalitz plot. The one- loop 
calculation does account for this phenomenon, but only to leading order in 
the low energy expansion. The final state interaction is analysed more accu- 
rately in two recent papers |]23| , p4| , which exploit the fact that analyticity 
and unitarity determine the amplitude up to a few subtraction constants. 
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Tl-» 71+71-71° (K + -K°) e , 

Q 




Figure 1: The l.h.s. indicates the values of Q corresponding to the various 
experimental results for the rate of the decay ?7^7r + 7r~7r . The r.h.s. shows 
the results for Q obtained with four different theoretical estimates for the 
electromagnetic self energies of the kaons. 



For these, the corrections to the current algebra predictions are small, be- 
cause they are barely affected by the final state interaction. Although the 
dispersive framework used in the two papers differs, the results are nearly 
the same: while Kambor, Wiesendanger and Wyler obtain r = 209 ± 20 eV, 
we get To = 219 ± 22 eV. This shows that the theoretical uncertainties of the 
dispersive calculation are small. 

Unfortunately, the experimental situation is not clear pSfl . The value of 
T^k+k-tt relies on the rate of the decay into two photons. The two different 
methods of measuring T r; ^ 77 (photon-photon collisions and Primakoff effect) 
yield conflicting results. While the data based on the Primakoff effect are 
in perfect agreement with the number Q = 24.2, which follows from the 
Dashen theorem, the 77 data yield a significantly lower result (see l.h.s. 
of fig. 2). The statistics is dominated by the 77 data. Using the overall 
fit of the Particle Data Group, r r) ^ 7r + 7r - 7r o = 283 ± 28 eV [p5 |, and adding 



errors quadratically, we obtain Q = 22.7 ± 0.8, to be compared with the 
result Q = 22 A ± 0.9 given in ref. |23|]. With this value of Q, the low 
energy theorem (|5|) implies that the electromagnetic self energy amounts to 
(Mk+ —MKo) e .m. = 2 MeV, to within an uncertainty of the order of 20%, in 
agreement with the lattice result. I conclude that, within the remarkably 
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small errors of the individual determinations, the two different methods of 
measuring Q are consistent with each other, but repeat that one of these relies 
on the lifetime of the 77, where the experimental situation is not satisfactory. 
8. Kaplan and Manohar [17j pointed out that a change in the quark 



masses of the form m' u = m u + am^m,, (cycl. u — > d — > s — > u) may be 
absorbed in a change of the effective coupling constants L 6 , L 7 , L 8 . The re- 
sults obtained with the effective Lagrangian for the meson masses, scattering 
amplitudes and matrix elements of the vector and axial currents are invari- 
ant under the operation. Conversely, since the ratios m u jm d and m s jm d 
do not remain invariant, they cannot be determined with the experimental 
low energy information concerning these observablesQ. In particular, phe- 
nomenology by itself does not exclude the value m u = 0, widely discussed in 
the literature p6j , as a possible solution of the strong CP puzzle. 



We are not dealing with a symmetry of QCD, nor is the effective La- 
grangian intrinsically ambiguous: even at the level of the effective theory, 
the predictions for the matrix elements of the scalar and pseudoscalar opera- 
tors are not invariant under the above transformation. Since an experimental 
probe sensitive to these is not available, however, the size of the correction 
Am in eq. (0) cannot be determined on purely phenomenological grounds - 
theoretical input is needed for this purpose. In the following, I use the 1/N C 
expansion and the requirement that SU(3) represents a decent approximate 
symmetry. For a more detailed discussion of the issue, I refer to ||. 

9. The problem disappears in the large- N c limit, because the transfor- 
mation m' u = m u + am d m s violates the Zweig rule || |[. For N c -^oo, the 
quark loop graph that gives rise to the U(l) anomaly is suppressed, so that 
QCD acquires an additional symmetry, whose spontaneous breakdown gives 
rise to a ninth Goldstone boson, the if. The implications for the effective 
Lagrangian are extensively discussed in the literature, and the leading terms 



in the expansion in powers of 1/N C were worked out long ago [27]. More 
recently, the analysis was extended to first non-leading order, accounting for 
all terms which are suppressed either by one power of 1/N C or by one power 
of the quark mass matrix ]2"B|| . This framework leads to the bound 



Am > , (8) 



lr The transformation maps the elliptic constraint onto itself: to first order in isospin 
breaking, the quantitiy 1/Q 2 may equivalently be written as (m 2 , — mfy/im 2 — ^d)> & n d 
the differences m 2 d — m„ , m 2 — m 2 , , — ml are invariant. 
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Figure 2: Quark mass ratios. The dot corresponds to Weinberg's values, 
while the cross represents the estimates given in ref. |§. The hatched region 
is excluded by the bound Am > 0. The error ellipse shown is characterized 
by the constraints Q = 22.7 ± 0.8, Am > 0, R < 44, which are indicated by 
dashed lines. 

which excludes the hatched region in fig. 2. Since the Weinberg ratios 
correspond to Am = 0, they are located at the boundary of this region. In 
view of the elliptic constraint, the bound in particular implies m u /md^ J \ 
and thus excludes a massless w-quark. 

10. An upper limit for m u /md may be obtained from the branching ratio 
T^r-yfap/T^i-^. Disregarding electromagnetic contributions |j29| , the ratio 
of transition amplitudes is proportional to {m d — m u )/{m s — rh): 

T^^^o 3y/3 , . 1 _ m d - m u 

l^'-^+fj 4 it R m s — m 

SU(3) predicts that, for quarks of equal mass, A^ vanishes: this term repre- 
sents an SU(3)-breaking effect of order m s — rh. The data on the branching 
ratio imply R = (31 ±4) (1+A^), where the given error only accounts for the 
experimental accuracy. The breaking of SU(3) is analysed in ref. ||29fl , on the 
basis of the multipole expansion. The calculation yields a remarkably small 
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result for Ay,/, indicating a value of R close to 31, but the validity of the 
multipole expansion for the relevant transition matrix elements is doubtful 
Moreover, fig. 2 shows that the result of this calcuation is in conflict 
with the large- N c bound. Since the quark mass ratios given in refs. H rely 
on the value of R obtained in this way, they face the same objections. 

At the present level of theoretical understanding, the magnitude of Ay/ is 
too uncertain to allow a determination of R, but I do not see any reason 
to doubt that SU(3) represents a decent approximate symmetry also for 
charmonium. The scale of first order SU(3) breaking effects such as Am, 
(F K - F w )/F w or Ay,/ is set by {M 2 K - M^)/M 2 ~ 0.25. Indeed, a correction 
of this size would remove the discrepancy with the large- N c bound. Large 
values of R, on the other hand, are inconsistent with the eightfold way. 
As a conservative upper limit for the breaking of SU(3), I use |Ay,/| < 0.4. 
Expressed in terms of R, this implies R < 44. The value m s /m = 29 ± 7, 
obtained by Bijnens, Prades and de Rafael with QCD sum rules [§], yields an 
independent check: the lower end of this interval corresponds to Am < 0.17. 
Figure 2 shows that this constraint also restricts the allowed region to the 
right and is only slightly weaker than the above condition on R. 

11. The net result for the quark mass ratios is indicated by the shaded 
error ellipse in fig. 2, which is defined by the following three constraints: 
(i) On the upper and lower sides, the ellipse is bounded by the two dashed 
lines that correspond to Q = 22.7 ± 0.8. (ii) To the left, it touches the 
hatched region, excluded by the large- N c bound, (iii) On the right, I use the 
upper limit R < 44, which follows from the observed value of the branching 
ratio ry/^y^o/Ty/^y,^. The corresponding range of the various parameters of 
interest is 

777 777 777 

— = 0.553 ±0.043 , — = 18.9 ±0.8 , — = 34.4 ± 3.7 , 

m d m d m u 

fin — nfh jYi 

40.8 ± 3.2 , = 24.4 ± 1.5 , A M = 0.065 ± 0.065 . 



m d -m u m 

While the central value for m u /m d happens to coincide with the leading 
order formula, the one for m s /m d turns out to be slightly smaller. The 
difference, which amounts to 6%, originates in the fact that the available 
data on the rj lifetime as well as the lattice result for the electromagnetic 
self energies of the kaons imply a somewhat smaller value of Q than what is 
predicted by the Dashen theorem, in agreement with ref. H. 
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The result for the ratio of isospin- to SU(3)-breaking mass differences, R = 
40.8±3.2, confirms the early determinations described in [Q]. As shown there, 
the mass splittings in the baryon octet yield three independent estimates of 
R, i.e. 51±10 (N-P), 43±4 (£"-£+) and 42±6 (~--H°)0. These numbers 
are perfectly consistent with the value given above. A recent reanalysis of 
p-Hjj mixing |3T| leads to R = 41 ±4 and thus corroborates the picture further. 



I find it remarkable that, despite the problems generated by the deter- 
minant of the Dirac operator for quark masses of realistic size, the lattice 
results for the mass ratios are quite close to the above numbers. The most 
recent values are m u /md = 0.512 ± 0.006, (m^ — m u )/m s = 0.0249 ± 0.0003, 
where the error only accounts for the statistical noise [[K|] . They correspond 



to Q = 22.9, Aju = 0, R = 38.6 - the place where the error ellipse shown in 
fig. 2 touches the large- N c bound. 

Finally, I use the value of m s obtained with QCD sum rules |7], [| as an 
input and calculate m u and with the above ratios. The result for the 
running masses in the MS scheme at scale fi = 1 GeV reads 

m u = 5.1 ±0.9MeV , m d = 9.3 ± 1.4 MeV , m s = 175 ± 25MeV . 
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